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Abstract 


Stationary spot welds have been made at the (001) surface of Fe-15%Ni-15%Cr single crystals 
using a Gas Tungsten Arc (GTA). On the top surface of the spot welds, very regular and 
concentric “ripples” were observed after solidification by differential interference color 
microscopy. Their height (typically 1-5 pm) and spacing (typically -60 pm) decreased with the 
radius of the pool. These ripples were successfully accounted for in terms of capillary-wave 
theory using the fundamental mode frequency fo given by the first zero of the zero-order Bessel 
function. The spacing d between the ripples was then equated to Vj/fo, where v s is the 
solidification rate. From the measured ripple spacing, the velocity of the pool was deduced as a 
function of the radius, and this velocity was in good agreement with the results of a heat-flow 
simulation. 
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1. Introduction 


The formation of ripples during welding or laser remelting is a problem that has been known for 
many years but which is yet not fully understood. Several mechanisms have been invoked to 
explain the ripple formation, most of them being related in some way to capillarity. Ripples are 
sometimes associated with what has been called “low-velocity bands” [1,2] : in this case, the 
whole melt pool does not seem to move with a steady-state speed and bands can be seen in 
transverse or longitudinal sections of the weids. In some other cases, ripples appear only at the 
top surface without in-depth correlation. 

Kotecki et al [3] studied the formation of ripples in GTA (Gas Tungsten Arc) spot welds. They 
observed concentric ripples after solidification which they attributed to oscillations of the liquid 
surface. However, their analytical model of natural oscillation in terms of a vibrating drum 
membrane is over-simplistic. More recently, Postacioglu et al [4] made a detailed calculation of 
natural and forced oscillations of a melt pool with particular attention to keyhole formation. 

In a now classical paper, Anthony and Cline [5] took a totally different approach to explain 
ripple formation during continuous laser welding. They considered the simple situation of a 
one-dimensional Marangoni flow induced by the surface tension gradient and calculated the 
pressure difference along the surface. They converted then this pressure gradient into a 
metallostatic head difference but ignoring the capillarity contribution (i.e., the Laplace force). 
Nakane [6] has made a detailed experimental investigation of ripples formed during GTA 
welds, under various conditions. He made a clear distinction between macroripples, which 
occur at low speed and propagate in depth of the weld, and microripples appearing at higher 
speed. 

In the present work [7], spot welds were made in austenitic single crystals in order to study the 
formation of microstructure and ripples. It is shown that the small ripples that form under these 
conditions can be analyzed using the formalism of standing capillary waves [8]. This formalism 
is summarized in section 3 after the experimental section. In section 4, the modeling of fusion 
and resolidification of the molten pool is briefly presented. Finally, the results are presented and 
discussed in section 5. • 


2. Experiment 

The spot welds were made on Fe-15 wt%Ni-15%Cr single crystal specimens. The procedure 
used to grow and to prepare these specimens is identical to that used to perform autogeneous 
welds [1] and thus will be only briefly summarized here. After being grown using a Czochralski 
technique, the single crystals were oriented by X-ray diffraction and then cut in order to obtain 
thin plates (typically 4 mm thick and 30 mm in diameter) bounded by (100) flat surfaces. These 
surfaces were then polished and ultrasonically cleaned in an acetone bath. 

The spot welds were made at the center of the discs using a Gas Tungsten Arc (GTA) process, 
under the protection of a He-Ar shielding gas. The standard conditions were typically 80 A and 
10 V applied for 10 s by a continuous power source. The resulting spot weld at the end of the 
heating stage was typically hemispherical in shape, with a diameter of about 6 mm. The power 
was then stopped and the molten pool solidified typically m about 0.5 s. 

After solidificauon, the top surface of the spot welds was directly observed without any 
preparation by optical microscopy, sometimes using a differenual interference color adapter. 
This last equipment could reveal very nicely the small ripples formed at the surface. The surface 
profiles were measured along radial lines using a profilometer with a diamond tip. 
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3. Capillary wave theory and ripple formation 


This section of the paper addresses the mathematical formulation of standing capillary waves 
and the resulting ripples that form during solidification. The model is derived from the 
formalism of Landau and Lifshitz [8] : it is based upon the Navier-Stokes equation with 
appropriate boundary conditions and is much simpler than the formalism derived by Postacioglu 
et al [4] using an energy principle. 

The basic assumptions of the model are as follows : i) the fluid is incompressible, non-viscous 
and irrotational (i.e., rot v = 0) ; li) the hemispherical melt pool of radius R and height h is 
replaced by a cylindrical shape with a flat bottom having the same dimensions (see Figure 1) ; 
iii) the wave is totally constrained at the bottom and lateral surface of the molten pool, 
including the triple point between the liquid, solid and air (i.e., zero velocity of the fluid at r = R 
and z = -h, see Figure 1) ; iv) only the vertical component of the velocity, v z , is considered. If 
these hypotheses are reasonable in the axisymmetric situation and the small amplitude wave 
considered here, the question of the constraint introduced at the triple point will be discussed 
after presenting the theory. Neglecting the viscous term, the Navier Stokes equation for the 
vertical component of the velocity, v f , can be written as : 

dv z 3v z 3v z 3P 

p 'dT=i r 3r + P v -'3;=-aI-p« <‘> 



Figure 1 : Schematic of the idealized 
situation used in the capillary wave model 
when the triple point is constrained not to 
move (a). The first mode which conserves 
mass under this condition corresponds to 
the second zero of the first-order Bessel 
function, J,,. If the triple point is allowed to 
move laterally (b), the first mode then 
corresponds to the first zero of J 0 


In this equation, p is the fluid density, P is the 
pressure and -g is the vertical component of 
gravity (assumed to be opposite to the z- 
direction). Since the fluid is irrotational (i.e., 
rot v = 0), this means that there is a velocity 
potential, <j), such that v = grad <{). Inserting 
this relation in eq. (1) and integrating with 
respect to z, one gets : 

dp 

P-P a = -pgz-p-^- (2) 


On the other hand, the difference between the 
dynamic pressure in the liquid, P, and the 
atmospheric pressure, P a , must be 
compensated by the surface tension of the 
liquid surface. Under the assumption of small 
displacements, this means that : 


P-Pa = 



1 3Z 
+ r 3r_ 


(3) 


where/ is the surface tension of the liquid 
and Z(r,t) the position of its top surface as a 
function of the radial distance and time. 
Please note that the curvature has been 
approximated by the Laplacian of the liquid 
position, -AZ. This equation can be derived 
with respect to time so as to introduce the 
vertical velocity of the top surface, and thus 



the velocity potential. Combining this equation with equation (2) also derived with respect to 
time, one finally gets an equation governing the potential <(>(r,z,t): 




p g 37 + 


3 2 <j> 
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at z = Z 


(4) 


Note that this equation is only valid at the top surface of the liquid and in reality is a boundary 
condition for the velocity potential <J>. The other equations governing the velocity potential are : 


1 d 3d>"| 3 2 <J> 

In the volume: . A0 = = 0 


(5) 


On the constrained surfaces of the cylinder, the velocity must be zero when a no-slip condition 
is imposed. Because only the vertical component of the velocity field, v z , has been considered, 
the following boundary conditions are imposed at the solid-liquid interface : 

3<j> 

v z = ~ = 0 at r = R and z = -h (6) 


A solution of Eqs (4-6) is provided in Ref. [7]. for the case of a cylindrical molten pool. In the 
case where the gravity term can be neglected in Eq. (4), the eigenmodes of the cylindrical 
molten pool are given by : 


(7) 

where Xo, n is the n-tK zero of the zero-order Bessel function, J 0 , and f n the associated oscillation 
frequency of the molten pool. The first mode of vibration (A^j = 2.405), called here the 
“breathing” mode, would correspond to a uniform displacement of the pool surface (i.e., all the 
points of the free surface go up and down at the same time). Therefore, this mode does not 
satisfy mass conservation if the triple point is totally constrained (i.e., if the point at the junction 
of the solid, liquid and air is not allowed to move slightly). However, when the wetting angle at 
the free surface of the pool is small, which is the case for most metals as the liquid wets the 
solid very well, lateral displacement of this triple point can be envisaged (see Figure 2b). 
Although the present model does not strictly apply to such a situation, it will be assumed that 
Eq. (7) still gives a reasonable estimate of the fundamental oscillation frequency of the molten 
pool. This point is further discussed in Ref. [7]. 

The boundary condition at the triple point being relaxed, this point can move back and forth 
(and up and down). In this case, the situation is similar to the basic oscillation mode of a wave 
that would break down near the “shore” of the pool, while at the same time the pool solidifies 
from the edge. This will finally leave ripples on the frozen surface separated by a spacing, d, 
given by : 

d(R)= lJ(i) = v s( R ) T °( R ) (8) 

where v s is the solidification rate. The basic mode of vibration, f 0 , given by the first zero of J c , 
has been assumed in the present equation, x 0 being the time period of the mode. This situation is 
illustrated in Figure 2. If the liquid surface had a steady shape, the triple point would move 
almost linearly. However, such is not the case due to the capillary oscillations. At time t, the 
liquid surface is assumed to have the maximum lateral extent, and thus the minimal angle, Q mn , 
with the surface of the solid. One period. T 0 , of the oscillation later, the triple point has moved a 



Table 1 Thermophysical properties of the Fe-15%Ni-15%Cr 


Volumetric latent heat of fusion 


1.9 10 9 Jm‘ 3 

Volumetric specific heat at 

T = 200 °C 

3.16 10 6 Jm‘ 3 


T = 1400 °C 

5.45 10 6 

Thermal conductivity at 

T = 200 °C 

14 Wrri 


T = 1405 °C 

28 


T > 1425 °C (liquid) 

280 

Liquidus temperature 


1425 °C 

Solidification interval 


20 °C 


5. Results and discussion 

Figure 3 shows the interference color picture of a weld made under standard conditions on the 
(001) surface of a single crystal stainless steel specimen. The diameter of the molten pool at the 
end of heating (10 s) was about 6 mm. This value is in close agreement with the result of the 
simulation which predicts a molten pool of 6.6 mm in diameter after heating. 

During solidification, nearly 
perfectly circular and concentric 
ripples form and appear as 
successive white and black rings 
in this macrograph. They 
correspond to slightly raised and 
depressed parts of the solidified 
surface, respectively. As can be 
seen, their spacing decreases as 
the radius of the pool decreases. 
Equation (8) predicts that the 
spacing of the ripples should vary 
as the product of the 
solidification rate times the 
period of the capillary wave. As 
will be shown in more detail later 
in this section, the solidification 
rate is fairly constant up to mid- 
distance of the pool and thus this 
dependence of the ripple spacing 
Figure 3 : Differential interference color picture of the is mainly due to the R 3/2 - 

top surface of a spot weld made under standard conditions. dependence appearing in Eq. (7) 

(providing the depth over width 
ratio of the pool, h/R, is fairly constant). Near the center of the weld, the ripples are so narrowly 
spaced that they can hardly be seen. Please note that the ripples are more regular in the bottom 
part of the pool, whereas they look slightly blurred in the upper one. Precisely in this last region, 
the edge of the weld does not seem very well defined : this is due probably to the reminiscence 
of the heating stage which clearly induced some fluctuations in the weld pool position 
(instability of the arc -pool system). 

From the contrast seen in Figure 3, it also appears that the overall shape of the solidified surface 
is significantly higher near the periphery of the spot weld, then goes down close to the center 
(dark wide nng located near R = 1 mm), then goes up again (light ring). A cavity seems to have 
formed at the center (dark circle). This overall variation of the solidified surface has been 
analyzed in Ref. [7] in terms of solidification shrinkage. Superimposed to the contrast 




distance d = v s t 0 . At the half penod, to/2, the pool has the minimal lateral extent (or maximal 
pool height at the center) and the angle between the liquid and solid surfaces, &max, has now a 
maximum value. As a consequence, the triple point is slightly depressed with respect to the 
position it would have (dashed line) if the liquid surface shape was stationary. The trajectory of 
the triple point has a sinusoidal shape instead of a straight line as indicated in this figure. In 
summary, the solidified surface near the nm of the pool has a maximum height when the liquid 
pool has a minimum level at the center, and vice versa. 


Since both v s and f 0 are a function of the radius of the pool, R, the spacing between the npples 
also depends on R. As can be seen in Eq. (8), it is necessary to deduce a relationship between 
the growth rate and the radius of the pool. The next section outlines a heat-flow calculation 
based upon a purely-diffusive model with latent heat release which allows to obtain this 
velocity. 



Figure 2 : Schematics of the mechanism involved in the formation of ripples, 
with the vibration mode of Figure lb. 


4. Hemispherical pool fusion and solidification 

The formation of the molten pool and its subsequent resolidification was modeled using the 
finite element solidification code 3-mos, developed at the Ecole Polytechnique Federale de 
Lausanne [9]. Although this package can include convection, solute transport or microstructure 
calculations, only the heat-flow part was used. The heat flow equation with latent heat release 
was solved for an axisymmetric situation using an enthalpy method. 

Since the small specimen was attached onto a base plate of copper to ensure good electrical and 
thermal contact, the simulation was instead performed for a semi-infinite domain (in the present 
case a cylinder 20 cm in diameter and 4 cm thick). This avoided the problem of specifying an 
arbitrary heat transfer coefficient between the plate and the copper chill. The power (800 W) 
was distributed evenly over an arbitrary surface of 4 mm of diameter and applied for a period of 
10 s (the experimental heating period). The thermophysical properties of the specimen used in 
the heat flow computations are indicated in Table 1. 

The computation was performed by setting adiabatic conditions at the external boundaries, 
except of course under the arc. As the diffusion distance over a period of 10 s is only about 
7 mm, setting a non-zero heat transfer coefficient on these surfaces has no effect. Similarly, 
radiation and convection on the welded surface are negligible with respect to the heat sink 
provided by the copper chill (or the bulk of the specimen in the present case). The initial 
temperature was set to room temperature. The velocity at the top surface was then calculated 
from the evolution of the molten pool shape 


associated with the uneven surface of the weld are faint straight lines at ±45 deg with respect to 
the horizontal. They correspond to the [±100] and [G±10] directions of the dendrite trunks 
growing nearly parallel to the surface. As was shown for autogeneous welds made on single 
crystals [1], the dendrite trunks which are selected by the growth competition process are those 
which have the minimal velocity or minimal undercooling. In the case of the non-stationary spot 
welds investigated here, the thermal gradient and velocity of the liquidus isotherm are 
approximately radial. Therefore, within each quadrant of the solid-liquid interface defined by 
the (110) and (110) planes, a specific <100> dendrite trunk variant grows near the top surface. 

The profilometry measurement of the weld shown in Figure 3 is presented in Figure 4. It has 
been recorded for a region near the periphery of the weld, along a radial direction, and the scales 
are indicated in the figure caption. As already observed on the macrograph of Figure 3, the 
spacing between the ripples becomes smaller as the radius of the pool decreases. The amplitude 
of the ripples is typically 1 pm and also decreases with the radius. The overall shape of the 
surface is remarkably flat, less than 5 pm variation over the 1.6 mm analyzed distance. 



Figure 4 : Profilometry measurement along a radius of the weld shown in 
Figure 3. The horizontal scale is 1.6 mm starting from the edge of the 
weld, the vertical displacement measured on the surface is 5.6 pm. 

Using the ripple spacing shown in Figure 4 and the model of capillary wave presented in section 
3 (Eq. (7)), the actual velocity of the pool can be computed as a function of its radius (Eq. (8)). 
The breathing mode (A*,, i = 2.405) was used for a “hemispherical” pool (h/R = 1) with a surface 
tension y= 1-9 N/m. The results obtained in this way are shown in Figure 5 (filled squares), 
together with another set of measurements (filled triangles) performed on the same weld but 
along another radial direction, perpendicular to that used in Figure 4. As can be seen, the data 
are fairly consistent. A least-squares fit of the type v(R) = A R‘ n was then made independently 
for each set of data. The result is displayed with dashed lines in Figure 5 with an indication of 
the fitted parameters. The prediction of the heat flow calculation is shown with open circles in 
this figure. The computation which was done independently from the ripples analysis shows 
that the results are remarkably consistent. In both analyses, the several parameters were not 
adjusted to obtain this agreement, but were taken from the available literature. The continuous 
curve shown in Figure 5 is also a fit v(R) = A R' n made on the data points obtained from the 
heat flow calculation over the same range of radius. 

Near the center of the weld, the solidification rate calculated with the heat flow model does not 
follow a R‘ n dependence, but unfortunately the ripples are no longer visible in this region. 
Nevertheless, the solidification times, t s , calculated from the npples analysis (t s = 0.61s or 
0.54 s for series #1 and #2, respectively) or from the heat flow computation (ts = 0.59 s) are in 
very good agreement. In the first analysis (ripples), this time was obtained from the integration 
of the fit v(R) = dR/dt = A R' n , from 0 to 3 mm, whereas that mentioned for the heat flow 
analysis directly comes from the computation. 
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Figures : Solidification rate as a 
function of the radius of the molten 
pool, as deduced from the measured 
ripple spacing and the model 
described in section 3 The squares 
and diamonds correspond to two 
sets of measurements on which 
dotted curves of the type v = A R " 
have been fitted using a least- 
squares method. The continuous 
curve shown in this figure has been 
calculated independently with the 
simple model of weld pool 
solidification outlined in section 4 


Radius [mm] 


7. Conclusion 

Based upon the capillary wave theory, a simple analysis of ripple formation has been 
performed. It has been shown that the ripple spacing is the result of the product of the 
solidification rate and of the period of the melt pool oscillation The solidification rate deduced 
from the analysis of the ripple spacing has been found to be m good agreement with that 
deduced independently from heat flow computations. 
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